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ABSTRACT 

Recent high-resolution spectra of the Type la SN 2006X have revealed the presence of time- variable 
and blucshiftcd NalD features, interpreted by Patat et al. as originating in circumstellar material 
within the progenitor system. The variation seen in SN 2006X induces relatively large changes in the 
total NalD equivalent width (AEW ~ 0.5 A in just over two weeks), that would be detectable at 
lower resolutions. We have used a large data set comprising 2400 low- resolution spectra of 450 Type 
la supernovae (SNe la) obtained by the CfA Supernova Program to search for variable Na I D features. 
Out of the 31 SNe la (including SN 2006X) in which we could have detected similar EW variations, 
only one other (SN 1999cl) shows variable NalD features, with an even larger change over a similar 
~ 10-day timescale (AEW = 1.66 ± 0.21 A). Interestingly, both SN 1999cl and SN 2006X are the two 
most highly-reddened objects in our sample, raising the possibility that the variability is connected 
to dusty environments. 

Subject headings: circumstellar matter — supernovae: general — dust, extinction 
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1. INTRODUCTION 

Type la supernovae (SNe la) are thought to result from 
the thermonucle ar runaway in a carb on-oxygen white 
dwarf (WD) star (|Hovle fe Fowlerll960l ). having gro wn m 
mass to the Chandrasekhar critical mass ddiandra sekharl 
1 193 If) . How the WD actua lly approaches this critical 
mass is still debated (e.g ., iParthasarathv et "all 120071 : 
ITutukov fc Fedoroval I2007D . but viable hypotheses in- 
clude the accretion of material from a non-degenerate 
binary companion ( the "single-degenerate" scenario; 
IWhelan fc M[il)73h . or a merger process involving an- 
other white dwarf star (the "double-dege nerate" sce- 
nario; llben fc Tutin^[l^8llWebbinkl[l98l . 

While there is no direct empirical evidence for either 
scenario, observations suggest the existence of multiple 
progenitor channels. For instance, the distribution of 
SN la peak luminosities is different in early-type and late- 
type galaxies (the most lu minous SNe la are f ound only 
in lat e-type galaxies; e.g.. lHamuv et al.1l2000t [Tha et al.l 
2007). Several studies have modeled the SN la rate 
with two components: one proportional to the star for- 
mation r ate, the other proportional to the total stel- 
lar mass dMannucci et al.l 120051 : IScannapieco fc BildsterJ 
120051: ISullivan et al.ll2006l ). The dependence on the star 
formation rate is evidence that at least some SNe la orig- 
inate in young stellar populations, and confirms earlier 
claims of an enhanced SN la rate i n star- forming galax- 
ies Ce.g.. lOemlerfc Tinslevi 119791: Ivan den BergH 119901: 
ICappellaro et al.l 119971 ). More recently. iPritchet et al.l 
( 2008) find that the SN la rate is ~ 1% of the WD for- 
mation rate, independent of star-formation history, and 
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conclude that the single-degenerate channel alone is not 
a viable model for all SN la progenitor systems (see also 
lYungelson fc Lividl200Ct lGreggidl2005h . 

The double-degenerate scenario has been invoked as 
a plausible interpretation for a fe w unusually lumi- 
nous SNe l a (e.g., SNLS-03D 3bb, iHowell et all 120061 
SN 2006gz, iHicken et ail l2007fh However, several hy- 
drodynamical models suggest that double- WD mergers 
cause an off-center carbon ignition that l eads to the con- 
versio n of C+O into O+Ne+Mg (e.g., ISaio fc Nomotol 
l2004h . and subsequently to an accretion-induced colla pse 
to a neutr on star (e.g.. iBravo fc Garcia-Send [l999t al- 
though see lYoon et al. I2007D . Radiation-hydrodynamics 
simulations of such O+Ne+Mg cores predict a weak 
<10 50 erg) explosion if rotation is ab sent, powered 
primarily by neutrino-en ergy deposition (|Kitaura et alJ 
l2006t [Dessart et "aLll2006|) . However, angular momentum 
is expected to be abundant in such accretion-induced col- 
lapses, resulting in a strong asphcrical explosion. The 
predicted low 56 Ni yields and ejecta mass would make 
these events op tically dim, quite unlike SN la explosions 
([Dessart et al.ll2007[ ). 

A direct confirmation of the single-degenerate sce- 
nario would be the detection of circumstellar material 
(CSM) arising from the transfer of matter to the WD 
by its non-degenerate binary companion. The interac- 
tion of the SN la ejecta with the CSM would proba- 
bly lead to detectable radio and X-ray emission, but 
obs ervations at these wavelengths yield only upper lim- 
its dPanagiaiEan [Hi EugEsiLin [2003) . Nonethe- 
less, strong hydrogen em ission has been dete cted in two 
SNe la so far (SN 2002ic,lHamuv et all [20031 : SN 2005gj, 
Aldering et all 120061 : iPrieto et al.l 120071: although see 



Trundle et all 2008 for a different 



Benetti et al.ll2006l and 

interpretation), interpreted as evidence for ejecta-CSM 
interaction, a nd hence for a sing le-degenerate scenario 
(although see lLivio fc R icss 2003|). However, the lack of 
similar detections in other SNe la suggests that these are 
rare events. 

Moreover, models suggest that a significant amount 
of material is stripped from the companion's envelope 
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during the exp losion (IChugail 1 19861 iLivne et all Il992t 
iMarietta et~al"T 120001 : iMeng et all l2007f ). which should 
lead to an observable signature in the form of a narrow 
Ha emission line in the ne bular spectra. No such signa- 
ture has yet been detected (jMattila et al.ll2"o"05t iLeonardl 
2007). A recent hydrodynamical investigation of the in- 
teraction between the ejecta of a Ty pe la supernova and , 
a main sequence binary companion ([Pakmor et al.ll2008ft 
finds that the theoretical predictions for the mass of 
stripped hydrogen are consistent with the observational 
limits. 

In a recent study, iPatat et"aT1 (|2007bf ) reported on the 
possible detection of circumstellar material in the normal 
Type la SN 2006X. Using high-resolution spectra at four 
epochs out to ~ 2 months past explosion, they observed 
a complex evolution in blueshifted absorption features 
associated with the NalD doublet. Because of the lack 
of similar evolution in the CallH&K doublet, they con- 
cluded that the time evolution of the Na I D complex is 
not due to line-of-sight effects in interstellar gas. Instead, 
the variation was attributed to ionization changes in cir- 
cumstellar material (CSM) within the progenitor system, 
due to radiation from the supernova. The expansion ve- 
locities of this CSM (~ 50 km s _1 ) were compatible with 
wind velocities for early red-giant stars, thereby favoring 
a single-degenerate progenitor scenario for SN 2006X. 

More recent studies of SN 2000cx and SN 2007af us- 
ing similar multiepoch high-resolution data have not 
detected any variable NalD features. No NalD ab- 
sorption is detected in SN 2000cx to the level of a 
few mA (|Patat et al.l l2007ah . setting an upper limit of 
2 x 10 10 cm~ 2 on the column densi t y of t he absorbing 
material. In SN 2007af. lSimon et all |2007f ) place an up- 
per limit of ~ 10mA on the EW of circumstellar NalD 
lines, corresponding to a column density of 6 x 10 10 cm -2 . 

It is therefore still unclear whether the variability in 
the NalD line is a common fea ture am o ngst SNe la, 
or if SN 2006X is a unique case. IChugail 12 008) argues 
that the structure observed by IPatat et all ( 2007bl ) can- 
not originate in circumstellar material, based on the ex- 
pected low optical depth of NalD (r < 10~ 3 ) in the 
spherically-symmetric wind of a typical red giant (RG) 
star. 

With a large sample of multiepoch high-resolution 
spectra of SNe la, one could address these questions di- 
rectly. However, such data require dedicated programs 
on 8-10 m telescopes, and signal-to-noise ratio (S/N) 
constraints restrict this study to bright events (my < 15 
at peak). 

In this paper, we instead use a large database of 
low-resolution spectra obtained through the CfA Super- 
nova Program to investigate changes in the Na I D dou- 
blet with time. The variations in the NalD profile in 
SN 2006X induced changes in its equivalent width (EW) 
of > 0.5 A between —2 and +14 d from maximum light. 
We have a unique set of hundreds of SNe la whose spec- 
tra we can examine for EW variations like those seen in 
SN 2006X. The statistical power of such a large sample 
enables us to place meaningful constraints on the fre- 
quency (and nature) of variable Na I D features in SNe la. 

We show that SN 2006X is not the only supernova 
with variable Na I D absorption, and the nature of the 
variable absorption we detect leads us to think that the 



explanation given by IPatat et al.l (|2007bf ) for SN 2006X 
may not apply in all cases. What's more, the absence of 
variable absorption in most of the SNe la for which we 
have adequate data is an important clue to the nature of 
SN la progenitors. 

The paper is organized as follows. In §[2]we present the 
spectroscopic sample used in this study. We use a Monte 
Carlo simulation to estimate the systematic error in our 
EW measurements in § [3l The results are presented in 
§ 0] and discussed in § [51 Conclusions follow in § [51 

2. SPECTROSCOPIC DATA 

To investigate the variability of NalD features in 
SNe la we have used a large spectroscopic data set 
obtained through the CfA Supernova Program (P.I.: 
R. Kirshner). Since 1994, we have obtained ~ 2400 opti- 
cal spectra of ~ 450 low-redshift (z < 0.05) SNe la with 
the 1.5 m Tilli nghast telescope at FLW O using the FAST 
spectrograph ([Fabricant et al.l 11998). The setup used 
(300-line/mm grating and 3" slit) yields a typical FWHM 
resolution of ~ 6 A (~ 300 km s _1 at NalD). Several 
spectra were pu blished in studie s of specific supernovae 
(e.g., SN 1998bu: IJha et al.lll999l ). while 432 spectra of 32 
SNe l a have recently been published bv lMatheson et all 
( 2008). We also have comple mentary multi-band optical 
photometry for most obje cts (jRiess et al.f l999: J h a et al.l 
2006; Hic ken et al.ll2009h . All published data are avail- 
able via the CfA Supernova Archive 6 . 

In SN 2006X, the bulk of the variation in the NalD 
EW occurs before ^ 15 d past maxim um light, with a 
timescale of — 10 d (|Patat et al.ll2007bl ). In order to in- 
vestigate the occurrence of similar cases in our data, we 
only consider objects for which we have at least one spec- 
trum before +10 d past maximum light, and for which 
the spectroscopic coverage spans at least 10 d. The re- 
sulting sample comprises 1826 spectra of 168 SNe la. 

All the spectra were obtained with the same telescope 
and instrument, and r educed in a consistent manner (see 
iMatheson et al.|[2008l for details). The uniformity of this 
data set is unique and enables an accurate estimate of 
systematic errors in our measurements. 

3. MONTE CARLO ERROR ESTIMATION 

Detecting small Na I D equivalent-width variations is a 
challenging task with low-resolution spectra. At the res- 
olution of the FAST spectrograph (FWHM sa 6 A), both 
Di (5895.92 A) and D 2 (5889.95 A) lines comprising the 
NalD doublet are blended together. Any variation in 
the strength of (wea ker) blueshifted com ponents (as ob- 
served in SN 2006X; IPatat et al.ll2007bh will then result 
in a small fractional change in the total (D1+D2) EW of 
a single NalD absorption feature. 

In what follows we use a Monte Carlo simulation to 
determine the minimum detectable EW variation as a 
function of resolution and signal-to-noise ratio (S/N). We 
use the results of this simulation to estimate the num- 
ber of objects for which we would detect EW variations 
comparable to those seen in SN 2006X. 

3.1. Simulation presentation 

6 http:/ /www. cfa.harvard.edu/supernova/SNarchivc. html 
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TABLE 1 

Simulation parameters 



Parameter 



Range 



Doppler parameter ( km s _1 ) 
Optical thickness of main NaID2 line a 
Number of Na I D sub-components 
Optical thickness of sub-components 3 
Blucshifts of sub-components ( km s _1 ) 
SN age (days from _B-band maximum) 
Resolution 1 " (A) 
Signal-to-noise ratio (per pixel) 



6 = 10 
— 1 < log ro < 5 

1 < "sub < 5 
-2 < log T , sub < 

< fsub < 100 
-15 < t SN < +85 

0.1 < FWHM < 15.0 

1 < S/N < 200 



" to is the optical thickness at line center. b We assume there are 
4.2 pixels per FWHM. 



We generate synthetic Na I D profiles consisting of one 
main doublet and several blucshiftcd sub-components. 
In the case of SN 2006X these were interpreted as be- 
ing of inters t ellar a nd circumstellar origin, respectively 
(jPatat et al.l l2007bf ) . The simulation parameters are 
summarized in Table [JJ and explained below. 

We start by generating the main Na I D doublet at high 
resolution (FWHM = 0.1 A), assuming pure absorption 
and a Voigt profile shape. The width of the profile is set 
by the so-called Doppler parameter b, which is a mea- 
sure of the thermal and turbulent motions within inter- 
stellar clouds. In all that follows we set b = 10 km s _1 
(this choice has negligible impact on our results). We 
parametrize the doublet according to the optical thick- 
ness of the stronger D2 line, To- (The optical thickness 
of the Di line is simply tq/2 by virtue of the ratio of 
oscillator strengths). 

We then generate a random number of up to five sub- 
components in a similar fashion, blueshifted by up to 
100 km s _1 with respect to the main component. This 
corresponds to the velocity ran ge of winds of early re d 
giant or late subgian t stars (e.g., Judge fc Stencelll991f h 
that was invoked bv iPatat et al . (2007b) to explain the 
observed blueshifts of the various Na I D sub-components 
observed in SN 2006X. In principle, these components 
could form in faster winds (e.g., from main sequence 
or compact helium stars), or through the interaction of 
remnant shells of successive novae with the wind of the 
donor star, in which case the blueshifts could be up to 
an order of magnitude larger. However, we deliberately 
limit our investigation to match the only observed case 
of SN 2006X. 

The top panels of Fig. Q] shows example high-resolution 
synthetic line profiles. In all cases the stronger main 
NalD components are identical (with log To = 1), while 
the strength of the three (weaker) blucshiftcd compo- 
nents increases from left to right, resulting in an increase 
of the total EW from ~ 1.7 to ~ 2.1 A. Note that, in our 
simulation, r spans six decades (Table [lj , such that we 
also investigate cases where the main Na I D component 
is weaker than the blueshifted sub-components. We fur- 
ther consider cases where the sub-components are at the 
same wavelength as the main component ( "blueshift" of 
km s _1 ), such that the main NalD lines vary while 
the blueshifted sub-components do not. This situation is 
particularly relevant to an example of variable Na I D we 
present in § 14.11 

In measuring equivalent widths the definition of the 
underlying supernova pseudo-continuum is a potential 




5880 5890 5900 5880 5890 5900 5880 5890 5900 

Rest Wavelength [A] 

Fig. 1. — Top: Example input normalized line profiles, con- 
sisting of the main NalD doublet and three weaker blueshifted 
sub-components. The strength of these components increases from 
left to right. The number in the lower-left corner is the equivalent 
width (in A) of the entire profile. Middle: As above but at the 
resolution of our spectroscopic data set (FWHM Ri 6 A). Note the 
different ordinate range. Bottom: As above with added Poisson 
noise (S/N = 60 per pixel). We also report the equivalent width of 
the degraded profile (with a typical error ~ O.lA; see q 13.21 . 



source of systematic error. We account for this by mul- 
tiplying the n ormalized Na I D p rofiles by SN template 
spectra (from [Hsiao et ail [2(307\l at ages between —15 
and +85 d from maximum light. The age distribution 
is matched to that of our data set. Over this time inter- 
val the underlying pseudo-continuum at the Na I D posi- 
tion changes significantly. In § 13.21 we show that this has 
negligible impact on the measured equivalent width. 

The result is then degraded in resolution through con- 
volution with a Gaussian profile, and resampled to 4.2 
pixels per FWHM resolution to match our data. We al- 
low for resolutions as low as 15 A, but illustrate the typ- 
ical resolution of our spectroscopic data set (FWHMw 
6 A) in the middle panels of Fig. [TJ While we are un- 
able to fully resolve individual components of the Na I D 
doublet, the overall increase in the total EW is clearly 
visible. 

Last, we degrade the S/N by adding Poisson noise to 
the lower resolution spectrum. Again, the overall in- 
crease in the total NalD EW is visible in the bottom 
panels of Fig. [TJ 

We then measure the equivalent width as follows: we 
define continuum regions on either side of the D2 and Di 
lines (or of the single Na I D line when the doublet is un- 
resolved) spanning a maximum of 20 A. We assume the 
line to span three times the width of the Doppler (Gaus- 
sian) core, or 3 x FWHM/2%/21n2 « 1.3 x FWHM, plus 
an additional ~ 2 A on the blue side to account for the 
possible presence of blueshifted sub-components. The 
continuum is then determined by fitting a second-order 
polynomial to the continuum regions, and the equivalent 
width is computed over the line span defined in the usual 
way. 

3.2. Simulation results 
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Fig. 2. — EW measurement error as a function of S/N for three 
FWHM resolutions. The dash-dotted curve corresponds to the 
resolution of our data set (FWHM 6 A). The ordinate axis on the 
right indicates the corresponding 3cr error on the difference between 
EW measurements. The horizontal dashed line corresponds to a 
minimum detectable EW difference of 0.5 A. 

The total equivalent width measured on the synthetic 
profile degraded in resolution and S/N is compared with 
the known input total EW. We then investigate the vari- 
ation of the mean and rms of the difference as a function 
of S/N. 

We summarize the simulation results in Fig. [2] for the 
case where log To = 1. We show the rms residuals (or 
measurement error) of the difference between the true 
and measured total Na I D EW for three FWHM resolu- 
tions (0.2, 6.0, and 12.0 A) as a function of S/N. As ex- 
pected, the error decreases with increasing S /N and reso- 
lution. For a 3cr detection of an EW variation of > 0.5 A, 
as seen in S N 2006X between — 2 and +14 d from max- 
imum light (|Patat et al.ll2007bt ). the measurement error 
must be < 0.1 A, which in turn implies S/N > 50 per 
pixel at a FWHM resolution of 6 A. Applying this S/N 
cut to our data further reduces the sample presented in 
§ [2] to 294 spectra of 31 SNe la. 

We find that the strength of the main Na I D compo- 
nent has a small systematic effect on our measurement 
accuracy. The mean difference between the input and 
measured EW is < 0.04 A for logr < 4 and S/N > 10 
per pixel (at a resolution of FWHM = 6 A) , independent 
of the number and strength of the blueshifted NalD 
sub-components. The shape of the underlying pseudo- 
continuum at Na I D (which changes significantly with 
the age of the supernova) contributes no additional er- 
ror. The measurement error due to resolution alone (i.e. 
the rms residuals for FWHM = 6 A at infinite S/N) is 
~ 0.01 A. In what follows, we assume a total systematic 
error <r sys = 0.05 A. 

4. RESULTS 

Our final sample comprises SNe la that span a large 
range of known light-curve properties. In pa rticular, 
the li ght-curve decline rate parameter Arai5 (B) ((Phillips 
fl993h ranges from .83 for the overluminous SN 1995al 
(|Phillips et~aLlll999h to 1.48 for the slightly underlumi- 
nous SN 2005am (F olatelli fc Ham uy 2008, private com- 
munication; see also lLi et al.ll2006T ). However, none of the 
SNe la in our sample has spectra or light curves similar to 
the more underluminous SN 1991bg (Ami5(i?) = 1.93; 



iPhillips et aflH999l) . 

We have measured the total Na I D EW on all 294 spec- 
tra using the same technique as for our simulated data. 
For each measurement, we use the corresponding error 
spectrum to derive the associated statistical error, which 
is added in quadrature to the systematic error derived 
from our simulation (er sys = 0.05 A; see § 13. 2|) . 

We take special care to include only the host-galaxy 
Na I D line when measuring the total EW. For redshifts 
cz < 800 km s _1 , the Galactic NalD line contributes 
to the host-galaxy line region defined in § [3| This af- 
fects only one SN la in our final sample, SN 1994D, 
whose host galaxy NGC 4256 is at cz = +592 km s _1 
(|Falco et al.lll999l) . However, iKing et alj (jl995l) report 
a recession velocity of +880 km s _1 at the SN position, 
which is the value we assume in this study. Nonethe- 
less, the presence of Galactic high-velocity clouds in- 
ferred from high- r esolut i on spectra (up to +250 km s — x ; 
IHq fc Filippenkol H9951 IKing et alJ fl995h causes some 
overlap in our spectra with the Na I D absorption in NGC 
4256 in the SN line of sight, and our EW measurements 
for SN 1994D are systematically biased towards higher 
va lues (EW = 0.32+° jj A , cf. 0.087+ 0.002 A reported 
bv lHo fc Filippenko|[l995l for the NGC 4256 component; 
see Table©. All other 30 SNe la are at cz > 800 km s -1 , 
and we interpolate over any Galactic Na I D for the EW 
determination. 

4.1. Variable NalD in the highly-reddened SN 1999cl 

The single example of variable Na I D lines observed to 
date, SN 2006X, is too slender a foundation on which to 
construct a general model for the variation. Nevertheless, 
we can test the null hypothesis of a constant EW against 
our data for each supernova. Rejecting this hypothesis 
constitutes detecting variability for an individual super- 
nova. 

Using a % 2 test for a constant EW model, we find that 
two SNe la in our sample display variable Na I D fea- 
tures: SNe 1999cl and 2006X. While variable sodium 
features have already been detected in h igh-resolution 
spectra of SN 2006X (jPatat et al.ll2007bf ). their detec- 
tion in our data gives us confidence in our methods and 
their detection in SN 1999cl constitutes a new result. 
The x 2 P er degree of freedom (dof) for the constant EW 
model is 11.23 for SN 1999cl and 4.47 for SN 2006X 
(Table©. We show the EW variation for both objects 
in Fig. © All other x 2 /dof are < 1.3. Interestingly, 
both SNe 1999cl and 2006X are amongst the three most 
highly-reddened objects in our sample (host-galaxy red- 
dening E(B - U)host > lmag; see § I4.2[) . The other 
highly-reddened object, SN 2003cg, shows no such varia- 
tion in the NalD EW (x 2 /dof = 1.07 for a constant EW 
model; Table©. The EW variation for this supernova is 
also shown in Fig. [3l The implications of this result are 
discussed at length in § [5l 

To our knowledge there are no (multiepoch) high- 
resolution spectra of SN 1999cl to confirm our find- 
ings, but the agreement between ou r measu rements on 
SN 2006X and those of iPatat et all (|2007bl ) illustrates 
the accuracy of our measurements and Monte Carlo er- 
ror estimation (see Fig.©. We have further checked that 
variations in seeing (and thus in spectral resolution for 
an object that does not fill up the 3" slit of the FAST 
spectrograph) do not correlate with the observed Na I D 
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Fig. 3. — Left: Time-evolution of the equivalent width of the NalD doublet for the three most highly-reddened SNe la in our sample: 
SNe 2003cg (filled triangles), 1999cl (filled squares), and 2006X (filled circles). Both SNe 1999cl and 2006X exhibit a variable NalD EW, 
while that of SN 2003cg remains constant over time. The larger filled circles at —2d and +14 d correspo nd to EW measurem ents on 
high-resolution (FWHM ss 7 km s" 1 , or ~ 0.14 A, at NalD) VLT+UVES spectra of SN 2006X published bv lPatat et~aH (I2007bl1 . Right: 
Normalized NalD profiles for SN 2003cg (top), SN 1999cl (middle) and SN 2006X (bottom), plotted in rest- frame wavelength. The black 
and red lines correspond to the smallest and largest EW, respectively. Note the difference in ordinate range, decreasing from top to bottom. 
The vertical dotted lines indicate the wavelength positions of the individual Di and D2 lines. [See the electronic version of the Journal for 
a color version of this figure.] 
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Fig. 4. — One minus the difference between the normalized 
NalD profiles with the smallest and largest EW, respectively, for 
SN 1999cl (top) and SN 2006X (bottom). The red line in the lower 
panel corresponds to the variation for SN 2006X inferred by de- 
grading the resoluti on of the high-resolution spectra published by 
IPatat c t al. (2007b]). The vertical dotted lines indicate the wave- 
length positions of the individual Di and D2 lines. The horizontal 
line in the upper panel shows the size of one the spectral resolution 
(FWHM Si 6 A). [See the electronic version of the Journal for a 
color version of this figure.] 



variations. 

Note that there is a significant EW variation inferred 
from high-resolution spectra at later ages i n SN 2006X 
(AEW « 0.36 A between +14 and +61 d; IPatat et all 
l2007bl) . However, this variation is below our typical de- 



tection threshold of ~ 0.5 A, and the data from our obser- 
vations at these later ages (when the supernova is fainter) 
do not generally meet our S /N requirements (see § [3|) . 

The maximum difference in the total NalD EW is 
1.66 ±0.21 A for SN 1999cl (i.e. 8.1a different from zero) 
and 0.83 ± 0.18 A for SN 2006X (4.5a). Without high- 
resolution spectra of SN 1999cl, we cannot place use- 
ful constraints on the change in Na I D column density 
this would represent. We note that the EW increases 
with time, and that the timescale for this variation is 
~ 10 d, for both objects (although this is a lower limit 
for SN 1999cl due to lack of suitable data). The EW 
variation for SN 2006X between the extrema EW values 
(i.e. between +3d and +12 d past maximum) is con- 
sistent with a linear variation, with % 2 /dof < 1 for a 
slope of ~ 0.08 Ad" 1 . This is not the case for SN 1999cl, 
with x 2 /dof > 2 for a linear fit. It is difficult to esti- 
mate the significance of the structure seen around — 4d 
for SN 1999cl with the available data. 

We can examine more closely the NalD profiles in 
Fig. [3] to study the nature of the variation seen in 
SN 1999cl, and in particular whether it could be different 
from that seen in SN 2006X . The high-r e solutio n data 
for SN 2006X presented by IPatat et all (|2007bh show 
that the variation entirely stems from weaker compo- 
nents blueshifted by ~ 50 km s" 1 with respect to the 
main saturated NalD components. Despite the lower 
resolution of our data, the Na I D profile with the largest 
EW (plotted in red in Fig. [3]) does indeed appear skewed 
towards bluer wavelengths. This is more clearly visible 
in the lower panel of Fig. |4] where we plot one minus the 
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TABLE 2 

Reddening and NaID EW measurements 
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band optical light curves usiii£ 


; the MLCS2k2 code of| Jha 
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EW. The upper and lower limits correspond to the maximum deviations from the weighted 
mean; col. (4): maximum EW difference. The lo" error appears in between parentheses; 
col. (5): AEW divided by its la error; col. (6): x 2 P cr degree of freedom for a constant EW 
fit; col. (7); number of degrees of freedom (simply the number of data points minus one); 
col. (8): age (in days from B-band maximum light) of the earliest spectrum; col. (9); age 
range (in days) of the spectra used in the fit. 



difference between the NaID profiles with the smallest 
and largest EW (black line). Degrading the resolution 
of the high- r esoluti on spectra of SN 2006X presented by 
iPatat et al.l (|2007bf ) reveals a profile also skewed towards 
the blue (red line). 

In contrast, the variation observed in SN 1999cl seems 
to affect equally all Na I D components between ~ 5885 A 
and ~ 5898 A (Fig.|U top panel). The largest variation is 
observed blueward of the Nal Di line (around ~ 5886 A, 
corresponding to a blueshift of ~ 200 km s _1 ). However, 
given the low resolution of our data (FWHM « 6 A, or 
~ 4 pixels), it is difficult to gauge the significance of the 
variation that is not centered at the rest wavelengths of 
the Na I D lines. Still, if all Na I D components (includ- 
ing blueshifted ones, if present) vary by a similar amount, 
this would suggest the variation is associated in part with 
the same material where the bulk of the reddening oc- 
curs, unlike SN 2006X where the weaker variable NaID 



components are not associated with the material causing 
the large reddening towards this object (see also § I5.2[) . 

Both SNe 1999cl and 2006X are part of the subgroup 
of SNe la with high velocity gradients in the Si II 16355 
absorption (HVG; see lBenetti et al. 2005). A larger sam- 
ple is needed to confirm the suggestion by iSimon et al.l 
(2007) that the variability is preferentially associated 
with the HVG subgroup. We note that our sample in- 
cludes several other HVG SNe la (e.g., SN 2002bo) that 
do not display variable NaID lines. 

Our measurements also show that not all blueshifted 
NaID components in SNe la are variable: IPatat et al.l 
(|2007bh noted that a high-resolution spectrum of 
SN 1998es around maximum light had similar NaID line 
profiles as SN 2006X around the same age, but we see 
no evidence for variability (EW = 1.6^q'^qA; see Ta- 
bled]). IPatat et"aT1 (|2007bf l also noted similar profiles in 
the overluminous SN 1991T, but this object is not part 
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of the sample presented in this paper. 

4.2. The relation between NalD equivalent width and 
host-galaxy reddening 

With NalD equivalent-width measurements for 31 
SNe la we can independently investigate its relation 
with host-galaxy color excess E(B — U)host, as has 
already been done by several authors using super- 
nova spectra (iBarbon et alj|1990t [Richmond et al.lll994t 
iTuratto et al.ll2003l ). In particular, we can test whether 
the Na I D EW can be used to estimate the host-galaxy 
reddening for objects with limited photometric data. 

The results are displayed in Fig. [5] and included in Ta- 
ble El Here the EW values are reported as the weighted 
mean of all measurements for a given supernova, while 
the error bars indicate the maximum deviations from 
the mean (this explains the large error associated with 
SN 1999cl, due to intrinsic variations rather than ran- 
dom error). The color excess in the host galaxy was de- 
termined from fits to mul ti-band optical l ight curves us- 
ing the MLCS2k2 code of I Jha et all (|2007h . corrected for 
the Gala ctic reddening i n the line-of-sight using the dust 
maps of iSchleeel etafl (fl998h . Over one half of these 
light curves are part of a new sample that will be pub- 
lished soon by the CfA Supernova Group (jHicken et al.l 

While there is a general trend of higher color ex- 
cess corresponding to higher NalD EW (such that 
E(B - U) hos t > 0.5 mag for EW > 2 A), the scatter 
in E(B — V)host a t a given EW is considerable. At 
EW « 2 A for instance, E(B - V) host ranges from ~ 0.2 
to ~ 1.5 mag. Any empirical relation based on SN la 
data can at best provide an approximate lower limit on 
E(B — y)host- The po ints do appear to clu ster around 
the bivariate relation of lTuratto et al.l (|2003h . albeit with 
a larger spread. The origin of the scatter could be due to 
a (global) variation in the dust-to-gas ratio in the host- 
galaxy interstellar medium, associated with variations in 
gas m etallicity and dust composition (e.g.. iDraine et all 
2007). Effects local to the SN circumstellar environment 
in the form of light echoes have also been found to im- 
pact the reddening derived assuming a standard color 
evolution for SNe la (|Wandl2005t iPatat et alJl2006h . 

We note that the relation of color excess with EW is 
expected to be non-linear for a fixed dust-to-gas ratio 
be cause of line saturat i on eff ects at high EW, as observed 
by iMunari fc Zwitterj (|1997f ) in Galactic O and B stars. 
Linearity at high EW could hold if the extinction takes 
place in a number of clouds at different velocities, each 
producing unsaturated features. 

5. CIRCUMSTELLAR OR INTERSTELLAR ORIGIN 
FOR THE VARIABLE NalD LINES? 

5.1. Conflicting interpretations of the NalD 
variability 

In SN 2006X, the variability o f bluesh ifted NalD lines 
was interpreted by IPatat et alj (|2007bl ) as the result of 
changing ionization conditions in circumstellar material 
within the progenitor system due to UV radiation from 
the supernova. The various sub-components in Na I D 
then originate in the pre-SN variable wind of an early red 
giant companion, or fr om remnant shells e jected in suc- 
cessive nova explosions (|Patat et al.l l2007b) . In the latter 



case, an asymmetric wind from the companion is needed 
to decelerate the material ejected during nova outbursts 
to match the observed < 100 km s _1 blueshifts of the 
NalD components. In principle, the low fraction of ob- 
jects for which we detect variable NalD lines (2/31, 
i.e. ~ 6%) could be due to viewing angle effects (see 
IPatat et al]|2007a| ). 

In both SNe 1999cl and 2006X, the NalD EW in- 
creases with age, as expected from Na II — >I recombina- 
tion. The variation occurs on a similar ~ 10 d timescale, 
although this is a lower limit for SN 1999cl due to lack 
of data (see Fig. [3]). At first sight, the Nal D variability 
seems c onsistent with the interpretation of IPatat et al.l 
(|2007bh . that it is associated with CSM in the SN la 
progenitor system. 

In a recent paper, IChueail ((20081) argues that the 
blueshifted Na I D components in SN 2006X cannot orig- 
inate within CSM material from a red giant wind. The 
optical depth in NalD would render it undetectable 
(even in high-resolution data), while that in CSM lines of 
CallH&K would be sufficient to detect them. The ob- 
served lack of vari ability in the s e calciu m lines is central 
to the argument of lPatat et alj (|2007bO that the features 
observed in Na I D cannot be due to intervening clouds of 
interstellar material that are not associat ed with t he su - 
pernova progenitor system. In contrast, Chugai (2008) 
use the same observational fact about Call to set an 
upper limit on the (spherically symmetric) wind from 
the RG companion. Unfortunately, we cannot study the 
variability in Call H&K lines with our data due to lower 
S/N in the blue portion of our spectra. 

The variable Na I D l ines o bserved in SN 2006X are 
interpreted by iChugail (|2008l ) as being due to clusters 
of interstellar clouds with different calcium abundances. 
Whether or not this is a plausible explanation is beyond 
the scope of this paper ; but a direct consequence of this 
(noted by Chugai 2008) is that the variability in NalD 
should be preferentially associated with highly-reddened 
SNe la. This is certainly the case in our sample: two of 
the three SNe la with E(B — U)host > 1 mag exhibit vari- 
able NalD. If these time- variant features were randomly 
associated to all SNe la, independent of their reddening, 
the probability of finding such features only in two highly 
reddened la is very low (~ 0.6%). This fact alone is not 
proof that the variability is related to interstellar absorp- 
tion, but is suggestive enough to warrant further scrutiny 
once multiepoch high-resolution spectroscopy for a larger 
sample is available. 

5.2. The Reddening Towards SNe 1999cl and 2006X 

Both SNe 1999cl and 2006X were found to be red- 
dened by dust with unusual values of Ry, the ratio of 
total to selective extin ction (Ry ~ 1.6 for SN 1999cl, 
Krisciunas etal] 120061 ; R v « 1.5 for SN 2006X, 
Wang et al.ll2008bD . to be compared with Ry sa 3.1 for 
the a verage value in the Milky- Way (e.g.. ICardelli et all 
1989). Spectropolarimetry of SN 2006X has also shown 
that the dust producing the pola rization is differe nt from 
the typical Milky Way mixture (jPatat et al.ll2008f) . Stud- 
ies of larger samples also conclude that dust in the host 
galax ies of SNe la is characterized by low v alues of Ry 
(e.g. JConlev et al.|[2007HHicken et aHl2009T ). 

However, not all SNe la reddened by non-standard 
dust display variable NalD lines: SN 2003cg has a 
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color excess E(B — VQ host ~ 1-1 mag and Ry ~ 1.8 
(El ias-Rosa et alj|2006l ). yet its NalD equivalent width 
remains constant from one week before to two weeks 
past maximum light (mean EW » 5.7 ± 0.2 A; see Ta- 
ble [2|) ■ This suggests the variability is not systematically 
related to non-standard dust. If low values of Ry arise 
in circumstellar material aroun d SNe la (see I Wan gj |2005l : 
iPatat et~al . 2006; Goobar 2008f), this in turn implies that 
the Na I D variability does not necessarily result from ef- 
fe cts local to the SN e nvironment . 

IPatat et"ail (|2007b[ ) attribute the large reddening to- 
wards SN 2006X to a dense molecular cloud, caus- 
ing the strongly saturated Na I D lines and a num- 
ber of weaker features (CN, CH, Cal, diffuse inter- 
stell ar bands, or DIB) at the same heliocentric veloc- 
ity (jCox & Patadl2008f l. The column density inferred 
from the strongly saturated Nal D lines is logiV(Nal) w 
14.3, which is two orders of magnitude larger than 
that of the most intens e time- variant features (log TV w 
12; IPatat et"afl [20 07b). In this picture, the redden- 
ing produced by the material where these variable fea- 
tures arise, if any, is negligible. Additional evidence 
comes from the detection of a light echo associated wit h 
SN 2006X (TWang et al.ll2008aH Crotts fc Yourdonl l200l. 
ICrotts k, YourdonT l 20081 ) present compelling evidence 
that this echo is generated in a dust sheet ~ 26 pc in 
front of the supernova, presumably where the bulk of 
the reddening occurs, and hence is not associated with 
circumstellar material. 

In SN 1999cl, the variation appears to affect all NalD 
components, including some with blueshifts of up to 
~ 200 km s~ 1 (see Fig. [4]), although the presence of such 
blucshifted components is difficult to confirm given the 
low resolution of our spectra (1 pixel corresponds to 
~ 75 km s _1 at NalD). We also note that we detect 



the DIB at - 6284 A in our low-resolution spectra of 
SN 1999cl at the same velocity as the NalD absorp- 
tion (the spectra of S N 1999cl have been published by 
iMatheson et alj 120081 and are publicly available at the 
CfA Supernova Archive). However, this feature is too 
weak (EW < 0.5 A) for us to confidently detect any vari- 
ability. 

The fact that we detect variable Na I D features only 
in the two most highly-reddened SNe la in our sample 
does suggest that the variability is somehow connected 
to dusty environments, and is not necessarily due to cir- 
cumstellar absorption. 

5.3. Variable Interstellar Lines in Galactic Stars, 
GRBs, and the Type Iln SN 1998S 

Time-variable interstellar absorption lines have been 
detected in high -resolution spectra of Galactic stars (see 
ICrawfordl 120031 and references therein). However, the 
EW variation is typically < 0.01 A and occurs on the 
timescale of years to decades (cf. ~ 1 A over ~ 10 d for 
SNe 1999cl and 2006X). Another important difference 
is that the variation occu rs in several interstel lar lines 
(e.g., both Nal and Call; ICrawford et al.ll2002f) as op- 
posed to Na I D only (as observed by IPatat "eTalll2007bl 
in SN 2006X). In most cases, the variation is associated 
with interstellar shells and bubbles (e.g., 7 of the 13 vari- 
able interstellar sight lines compiled by ICrawfordl 120031 
are associated with the Vela SNR). The timescale of the 
variation for these stars implies the exist ence of structure 
in th e ISM on scales of 10-100 AU (e.g. jLauroesch et al.l 
12000ft . 

Time- variable interstellar lines h ave also be e n obse rved 
in other extragalactic objects. lHao et all (|2007f ) re- 
port on the detection of variable Fell and Mgll lines 
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in low-resolution s pectra of GRB 060206 (although see 
iThone et all 12008). which they interpret as stemming 
from insterstellar patches of Mgll clouds of size com- 
arable to the GRB beam size (~ 10 cm). iD'Elia et all 
2008) have revealed variable lines of Fell in high- 
resolution spectra of GRB 0803 19B, whose fine structure 
was excited through pumping by the GRB UV photons. 
The distance to the absorbing material is constrained 
to be ~ 18-34 kpc away. (Note that the UV flux of a 
GRB is orders of magnitude larger than that of a SN la, 
whose ionizing/exciting capabilities are rather limited in 
space). In both cases the variation occurs over very short 
timescales (a few hours). 

As noted bv lPatat et alj (|2007bl ). the only other super- 
nova with variable (and blueshifted ) Na I D lines is th e 
(core-collapse) Type Iln SN 1998S (|Bowen et al.ll2000D . 
In this object, however, the presence of narrow P-Cygni 
profiles of H and He at the same velocity as the vari- 
able NalD lines is interpreted as the signature of out- 
flows from the supergiant progenitor, henc e associating 
the variability with cha nges in the CSM dBowen et al.l 
120001: iFassia et alj|2001t see also iFransson et alJl2005F 

6. CONCLUSION 

We detect variable Na I D lines in low-resolution spec- 
tra of 2/31 (~ 6%) SNe la: SNe 1999cl and 2006X. 
The S/N of our spectroscopic data set is sufficiently high 
(S/N> 50 per pixel) that we can reliably detect changes 
of > 0.5 A in the NalD equivalent width. This is com- 
parable to the total EW variation seen in S N 2006X be- 
tween —2 and +14d past maximum light (|Patat et al.l 
l2007bh . We conclude that either variable NalD fea- 
tures are not a common property of SNe la, or that the 
level of the variation is less on average than observed in 
SN 2006X. 

The limitations of this study due to the low resolu- 
tion of our spectra are overcome in part by the larger 
statistical power of our sample. 

Only a small fraction (~ 6% in our sample) of SNe la 
display variable NalD lines. The two SNe la for which 
we detect variable Na I D (SNe 1999cl and 2006X) are the 
two most highly reddened objects in our sample, suggest- 



ing that the variability could be associated with inter- 
stellar absorption. However, we detect no such variabil- 
ity in the other highly-reddened SN 2003cg. We cannot 
exclude with our data the possibility that a larger frac- 
tion of objects display EW variations below our detection 
threshold of ~ 0.5 A. 

Whether the material causing the large Na I D EW 
variations is connected to the SN la progenitor system 
ought to become clear from multiepoch high-resolution 
spectroscopy for a larger sample. The detection of such 
variable NalD features in an unreddened SN la would 
reduce the probability that the variability results from 
a purely geometrical effect in clouds of interstellar ma- 
terial. If variable Na I D lines remain a rare property 
amongst SNe la, this would either indicate that the vari- 
ability is not associated with circumstellar material (and 
hence to the SN la progenitor system) , or that there is a 
preferred CSM geometry for detecting the variability. It 
is premature with the current available data to place firm 
constraints on the nature of Type la supernova progeni- 
tor systems based on the presence or absence of variable 
NalD lines alone. 
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